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Abstract 

We compare vector boson fusion and quark antiquark annihilation production of 
vector boson pairs at the LHC and include the effects of anomalous couplings. Results 
are given for confidence intervals for anomalous couplings at the LHC assuming that 
measurements will be in agreement with the standard model. We consider all couplings 
of the general triple vector boson vertex and their correlations. In addition we consider 
a gauge invariant dimension-six extension of the standard model. Analytical results 
for the cross sections for quark antiquark annihilation and vector boson fusion with 
anomalous couplings are given. 


1 Introduction 

In this note we study vector boson pair production with possible anomalous couplings in 
proton proton collisions at the LHC. The motivation to study these processes has been 
twofold: 


1. If the electroweak symmetry breaking is not realized by a light Higgs boson, the sym¬ 
metry breaking will manifest itself by some strong interactions among longitudinally 
polarized gauge bosons |I|, 0. In general, the amplitudes for longitudinal vector bo¬ 
son scattering are then very large at high energies. Several models to describe the 
strongly coupled symmetry breaking, in particular the standard model with a heavy 
Higgs boson and technicolor inspired models, have been discussed [0 i 0 0- If an 
amplitude has been calculated within a specific model, a method to connect this am¬ 
plitude to parton parton scattering processes has to be employed. The conventional 
method [|^, |], |^, [T^ was to use the effective vector boson approximation (EVBA) |1II| . 
The EVBA was originally used only for longitudinally polarized vector bosons. It 
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was however also applied to all intermediate helicity states 0 despite of the known 
problems with the EVBA for the transverse helicities ||T^ . 

Avoiding the use of the EVBA, the quark-quark scattering processes qiq 2 
were calculated exactly in lowest order of perturbation theory |^, In addition to 
the vector boson scattering diagrams, diagrams of bremsstrahlung type have to be 
evaluated in the exact calculation. 


2. On the other hand one may assume that the symmetry breaking is realized by a light 
Higgs boson. In this case the dominating processes for vector boson pair production 
are those of direct quark antiquark annihilation, also called Drell-Yan processes. The 
rates for these processes are sensitive to the values of the couplings of the electroweak 
vector bosons among each other [O. Drell-Yan production with anomalous (=non- 
standard) couplings has been studied in ||I^-|20[ 

m- 


0{as) corrections have been taken 


into account in (2B[| . The vector boson scattering processes were not considered 
in these works. The common argument to omit these processes was that they are 
O(a^) and therefore suppressed with respect to the Drell-Yan processes. However, a 
particular case in which these processes can give a signihcant contribution is near a 
Higgs boson resonance. In the study |2^ of the signal of a resonant Higgs boson both 
the Drell-Yan processes, including the 0{as) corrections |^|, and the exact matrix 
element for qiq 2 —>■ were taken into account. Also in [^], the processes 

qiq 2 —^ were included. These calculations however were only for standard 

vector boson self couplings and the rates for the two different production mechanisms 
have not been explicitly compared. 

In summary, in the strongly interacting scenario particular attention was paid to the vector 
boson scattering processes while the analyses of vector boson self couplings only considered 
the Drell-Yan processes. 

Later on, the effects of various SU (2) l x U{1)y gauge invariant effective interaction terms 
among the electroweak vector bosons were investigated and the vector boson scattering 
processes were considered [^, ^ together with the Drell-Yan processes. It was found 
that the Drell-Yan contribution and the one of vector boson scattering were of comparable 
magnitude. However, as in p, the vector boson scattering processes were calculated using 
the EVBA for all intermediate boson helicities. 

Recently |^, we showed that an improved version of the EVBA can increase the 
reliability of EVBA calculations. In particular, the improved EVBA could well reproduce 
the result of a complete perturbative calculation for a process which is dominated by the 
transverse intermediate helicities. 

In this article we carry out a comparison of Drell-Yan production and vector boson 
scattering using the improved EVBA and including the influence of anomalous couplings. 
This work is thus a supplement to the existing analyses |]^ ^ ^ ^ in which the Drell- 
Yan processes have been considered in more detail (0(as) corrections were included and 
more rehned kinematical cuts were applied), but vector boson scattering was not discussed. 
We will study the general parametrization [^],|^|-|^ of the triple gauge boson vertices in 
terms of seven parameters, allowing for C- and P-violation. In addition, we will study an 
SU{2)l X U{1)y gauge invariant dimension-six extension of the standard model. Our work 
extends the works |^, in that all three C- and P-invariant gauge invariant dimension- 
six operators [^-[^ which affect the vector boson self-interactions are discussed. We 


note that the three C- and P-invariant trilinear couplings which potentially contribute to 
the experimentally relevant [|^ process of W^Z Drell-Yan production can be equivalently 
expressed in terms of the parameters of the three-parameter gauge invariant model. The 
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Figure 1: Diagrammatic representations for the production of a vector boson pair V 3 V 4 in 
the collision of two hadrons hih 2 - a: via the quark antiquark annihilation mechanism, b: 
via the 0{a‘^) parton reaction qiq 2 


same is true for the two C- and P-invariant couplings which potentially contribute to the 
similarly relevant process of hF^y production. 

In Section ^ we compare vector boson fusion and Drell Yan production in the three- 
parameter gauge invariant model. In Section we present parameter hts for the anomalous 
couplings which can be obtained from future LHC measurements assuming that standard 
model predictions will actually be measured. We discuss the full set of anomalous couplings 
and also give the unitarity limits for the set of couplings which we use. We also consider 
again the three-parameter gauge invariant model. In Appendix ^ we give analytical formulas 
for the cross sections for qq' annihilation into IF^y and W^W~ pairs in terms of the 

seven anomalous couplings. In Appendix ^ we give formulas for vector boson scattering 
cross sections for the gauge invariant model. 


2 Comparison of Vector-Boson Fusion and Drell-Yan 
Production 


To illustrate our results we calculate the invariant mass distributions of the cross sections 
for vector boson pair production at the LHC (pp collisions at = 14 TeV). We com¬ 

pute both the contribution from Drell-Yan production and from the 0{a'^) parton reaction 
Q 1 Q 2 which contains vector boson scattering. The two contributions are shown 

diagrammatically in Figure |I|. Both contributions are calculated in the Born approximation 
and we use the improved EVBA |^, ^ to calculate the latter contribution. We discuss all 
possible pairs of produced electroweak vector-bosons, W~^W~, ZZ, Zy 

and yy. We hrst present the results for the standard model and then for non-zero anomalous 
couplings. 


2.1 Calculational Procedure 
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2.1.1 Drell-Yan Production 


In the usual quark-parton description, the lowest order contribution conies from the Drell- 
Yan processes shown in Figure |l] (a). Three generic Feynman diagrams can contribute to 
any of these processes in lowest order (Fig. §). They correspond to the exchange of vector 
boson(s) in the s-channel and the exchange of fermions in the t- and in the w-channel. 
Only the vector-boson exchange diagrams receive a contribution from the vector boson self- 
couplings. The standard model differential cross sections for qq' V3V4 have been hrst 
given in [Q . The results for arbitrary aw can be found in . For arbitrary vector boson 


self couplings, demanding only Lorentz-invariance, the differential cross sections as well as 
the expressions for the helicity amplitudes have been recently given for all processes in 
analytical form in [^. We choose to repeat the formulas for the differential cross sections 
in Appendix ^ in a form in which the high energy behavior is immediately transparent. We 
note that the 0(as)-corrections to the lowest order cross-sections can be huge. For 
production they can reach up to 70% of the lowest order contribution and for 
production ||^ they can be even larger. Only the Born cross section will be considered 
here. 

The formula for the invariant mass distribution of the cross section for DjIA-pair pro¬ 
duction via gg'-annihilation in the collision of two hadrons hih 2 is given by 

-{hih2 qq' V3V4,, Shh)\ cut 
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{qq' ^ V 3 V 4 ). 


( 1 ) 


This formula is valid if either no cuts or a rapidity or a pseudorapidity cut on both produced 
vector bosons is applied. A pseudorapidity cut, in contrast to a rapidity cut, always excludes 
events near the hadron beam direction. In (|^), ^yshh and are the invariant masses of 

the hadron pair and the vector boson pair, respectively, y is the rapidity of the gg'-pair in 
the hih 2 c.m.s and x = My^y^/Shh- The quantities denote the parton distributions in 
the hadrons and the quantities Q\ are the factorization scales. 9 is the angle between the 
quark q and the vector-boson V 3 in the center-of-mass system of the quarks. Applying no 
cuts, the limits of integration in (|ID are |/max = |ln(l/a:) and 2 :min(l/) = = 1- We 

choose here to apply a pseudorapidity cut rj on the produced vector bosons. This cut is 
equivalent to a minimum required angle 'dmin between the direction of momentum of any of 
the produced vector-bosons and the hadron beam direction. The cut rj is related to 'dmin by 
tanh(? 7 ) = cos'dmin- The integration limits with an r^-cut in the hih 2 c.m.s are given by 
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( 2 ) 


and one has to require that Zmin < ^max- The upper sign of =F in (ID is for Zmin, the lower 
sign for Zmax- In (11), the quantity (5 = tanh(|/) is the boost-parameter for a transformation 
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Figure 2 : The generic Feynman diagrams for a process qq' —>■ V3V4 in lowest order of 
perturbation theory. 


from the {hih2) c.m.s into the (V3V4) {={qq')) c.m.s. and = tan^'dmin- Further we dehned 
7^ = 1/(1 — / 3 ^). The quantities = \Jq^ + and F^4 = \J+ Ml are the energies of V3 
and V4 in the {qq') cms, while M3 and M4 are their masses, q is the magnitude of the three- 
momentum of V3 or V4 in the (V3V4) cms-system. The last argument of the min-function 
which dehnes ?/max in ( 0 ) only plays a role near the threshold. In deriving (H) we assumed 
that the quarks have no transverse momentum with respect to the hadrons, but no other 
kinematical approximations were made. 

For large energies of the produced vector-bosons, q^ S> max(M|,M|), the limits of 
integration (|^) take on the simplified forms 


~ mm - In (1/a:), 7 
^max(|/) - -Zrain{y) - tanh 


{y- li/D- 


( 3 ) 


In this limit, the rj-cnt is identical to a rapidity-cut Y of the same magnitude. We choose 
a cut of the magnitude rj = 1 . 5 , corresponding to a minimum angle of = 25 °. For the 
relevant process pp + X the highest sensitivity to anomalous couplings is achieved 

with a cut of about this magnitude . 


2.1.2 Vector Boson Fusion 


The 0 {a'^) partonic reaction which is shown in Figure 0 (b) contains the vector-boson scat¬ 
tering processes V1V2 —> V3V4 as subprocesses. Three types of Feynman diagrams contribute 
to a generic process V1V2 —>■ V3V4. They correspond to vector boson exchange, a direct inter¬ 
action among the four vector bosons and Higgs boson exchange. Using the Feynman rules 


for the GIDS model given in we wrote the amplitudes as functions of the scalar prod¬ 
ucts of the external momenta and of polarization vectors. We evaluated them numerically 
without making further approximations. In Appendix P we give analytical expressions for 
the cross sections for W^Z —> W^Z, lU^y —>• W^Z and 7 —>■ IU^7 in a high energy 
approximation. Expressions for the amplitudes of these and other vector boson scattering 
processes can be found in [TT, ED- 

We calculate the invariant mass distribution of the cross-section for hih2 — >■ V1V2 —V3V4 
in the improved EVBA according to [Q . The formulas which have been given there apply if 
a rapidity cut is used. The corresponding expressions for a pseudorapidity cut are obtained 
by replacing Z min ft/). z^a,x{y) and yma.^ in 0 by the expressions (H). We use the exact vector 
boson pair luminosities of if V1V2 consists of two massive vector bosons. If a photon is 


involved, the Approximation 2 of with the photon distribution function of [Q is used. 
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2.2 Results in the Standard Model 


Figs. and ^ show the invariant mass distributions for all vector boson pair production 
processes in the standard model. We separately show the contributions from the processes 
hiV 2 —>■ V 3 V 4 and those from qq' —> V 3 V 4 . The mass of the Higgs boson was chosen to be 
Mj{ = 300 GeV. There is little effect (less than 15% of change in the contribution from vector 
boson fusion) on the results for and hF^y production if the mass of the Higgs boson 

is varied in between = 0.1 TeV and Mh = 0.8 TeV. The other electroweak parameters 
were chosen as a = 1/128, Mz = 91.19 GeV and Mw = 80.33 GeV. We use the Higgs boson 
width T// for the dominant decay modes into W^W~ and ZZ, Th = Th^w+w- + Th^zz, 
where 

f - , 9 . , . 

In (^), xv = 4:My/Mjj and = sin^( 6 *w), where 6w is the weak mixing angle, n is the 
number of identical particles in the state VV. For the parton distribution functions we use 
the set MRS(R2) which includes the latest HERA data and uses as(M|) = 0.120 as 
input parameter, a value favored by the LEP 1 data. A contribution from top quarks is 
neglected. For the scales Qf appearing as arguments of the parton distribution functions 
we use the quark-quark sub-energy, Qf = Sq^ For the elements of the GKM matrix we 
take |I4dP = iWsP = 0.95, |I4sP = iWdP = 0.05 and consequently assume no mixing of 
the third flavor with the other two flavors. If no GKM mixing is included at all none of the 
differential cross sections changes by more than 1 %. 

Figs. 1^ to ^ clearly show that the contribution from vector boson scattering is always 
an order of magnitude smaller than the contribution from gg'-annihilation (also if the sum 
over all V 1 V 2 is taken). The contribution may therefore indeed be neglected. 

Fig. I shows the ratio of the cross sections for pp —>• fF^( 7 , Z) — (W^Z produc¬ 
tion = the sum of W^Z and W~Z production, hF^( 7 , Z) intermediate states = the sum of 
IF ^7 and W^Z intermediate states) and for pp —>• qq' —>• W^Z as a function of My^Vi- The 
ratio of the integrated cross sectionsQis 12% (15%) for a cut of 7 = 1.5 (7 = 2.5). 


^If Approximation 2 of is used for the vector-boson distribution functions, = xspp has to be 
chosen instead, where x is the first argument of f^{x, Q^). 

^For this numerical evaluation we used Mh = 0.1 TeV in order to be able to compare with results in the 
literature. We integrated the cross sections between 0.5 TeV < Mwz < 2 TeV. 
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PP ^ + X at Vs ^ 14 TeV 



0 0.5 1 1.5 2 0 0.5 1 1.5 2 


Mv3vjTeV] 

Figure 3: The cross sections for W~Z, W^''^ and hF “7 production as functions of the 

invariant mass My^Vi of the produced vector boson pair for pp-collisions at ^/s = 14 TeV. 
Separately shown are the contributions from gV&’^iiihilation and from vector-boson fusion 
processes. A rapidity cut of ?7 = 1.5 has been applied. 
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PP ^ + X at Vs ^ 14 TeV 



0 0.5 1 1.5 2 0 0.5 1 1.5 2 


Mv3vjTeV] 

Figure 4: The cross sections for W^W~, ZZ,W^W^ and W~W~ production as functions 
of the invariant mass My^v^ of th® produced vector boson pair for pp-collisions at y/s = 14 
TeV. Separately shown are the contributions from gV^^^iiihil^-tion and from vector-boson 
fusion processes. A pseudorapidity cut of ?7 = 1.5 has been applied. 
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Figure 5: The same as Fig. ^ but for Z 7 and 77 production. 


We note that a different value of this ratio is obtained if the EVBA in leading logarithmic 
approximation (LLA) is used instead. In Bi the cross sections for pp —> hF^( 7 , Z) —>• W^Z 
and for pp ^ qq' —> W^Z were calculated and the LLA EVBA was used. Calculating the 
ratio of these cross sections, we obtain 57% for Y = 1.5 and 64% for Y = 2.5 for the case of 
a light Higgs boson (59% (Y = 1.5) and 65% (Y = 2.5) for Mh = 1 TeV). Likewise, if we 
repeat the calculation of |^, ^ (we used rj = 1.5, Mh = 0.1 TeV and integrated the cross 
sections in the region 0.5 TeV < Mwz < 2 TeV), we obtain a value of 52% for the ratio. 
For more details we refer to 


45 


These values of the ratio are thus much higher than the values obtained with the im¬ 
proved EVBA. The latter values are however in agreement with values following from , in 
which a complete (lowest order) calculation of the processes qiq 2 q'iq 2 W~^Z was carried out 
instead of an EVBA. Computing the ratio of the cross sections for pp —> giq '2 ^ q[q 2 W~^Z 
and pp —>• qq' W'^Z given in ^ one obtains 17% (21%) for Mh = 0.1 TeV {Mh = 1 
TeV). In summary, the improved EVBA calculation and the complete calculation both yield 
a value for the ratio which is between 10% and about 20%, while calculations using the LLA 
EVBA yield a value which is larger by more than a factor of 3. 

We remark that for Mh = 300 GeV even the Higgs boson peak (which is present only 
in W~^W~ and ZZ production) stays below the rate of gg'-annihilation. We finally note 
that the like-charge pair production processes pp —W^IV^ -|- X cannot proceed via qq' 
annihilation and might thus allow to directly observe vector boson scattering. 
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Figure 6 : The ratio of the cross sections for pp —hF^(Z, 7 ) — W^Z and pp — >■ qq' —>■ W^Z 
as a function of the invariant mass My^v^ y/^pp ~ TeV. 


2.3 Parametrization of Anomalous Couplings 

The model we use for the anomalous couplings was described in (GIDS model). In 
this model, the most general SU{2)l x f/(l)y-symmetric interaction terms of dimension 
six are added to the Lagrangian of the standard model. We restrict ourselves to C- and 
P-conserving interactions which contain no higher derivatives and explicitly contain vector 
boson self-interactions. There are three of those interaction terms which are described by 
the parameters avc, and They are related to the usual parameters x^, y.y and 

Sz, xz, yz, which parametrize the C- and P-conserving interactions of the ^W^W~ and the 
ZW^W~ vertex, respectively, by 

oz — — 22i.gi — -, Xz — —) —-(aw# + as$) —- x^, 

Sw SyrCw Sw CW Cw 


Cw , Cw 

yz — —^z — —aw, 

Sw Sw 


x^ 


— — aw# + a^^, y^y — — aw- 


(5) 


In (I) we also included the relations to the parameters Agff, Akz, Xz and Ak^, X^ of 


The reduction from the hve parameter case of Sz^xz^yz, and y^ to the three parameter 
case is manifest through the relations xz = —{swlcw)x^ and yz = {cw/sw)y-y which are 
implicit in (^. The three parameter model dehned in (j^) has already been obtained [42 


m 


from the assumption of a custodial SU{2) symmetry. The relation between x^ and 
Xz in (1^), Xz = —{sw/cw)x^, is a consequence of the exclusion of intrinsic SU{2) violation, 
i.e., of SU{2) custodial symmetry. The relation between y^ and yz in follows from the 
requirement of SU{2)l x P(l)y symmetry in the quadrupole interactions. In addition to 
trilinear interactions the three-parameter dimension-six SU{2)l x P(l)y gauge invariant 
model describes interactions among four and more vector bosons. Also these interactions 
are already contained in an identical form [Q in the model described in |T^ , ^ . The only 
difference of the three-parameter model |T^, ^ and the SU{2)l x P(l)y invariant one 


lies in non-standard interactions of the Higgs boson. 

We note that there are no non-standard interactions among three neutral gauge bosons 
which would obey C- and P-symmetry, contain no higher derivatives and are compatible 
with electromagnetic gauge and Lorentz invariance 
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'‘The parameters are called ew, ew# and es# in [M 





























The Lagrangian of the GIDS model is an effective, unrenormalizable one and can in 
general be written as R 


-^efr — -Cq 




9j ^(6) 

Z. A2^i 
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+ ... 


( 6 ) 


In (P), Cq is the Lagrangian of the standard model, the are interaction terms of dimen¬ 
sion d, the Qj are conpling constants and A is the energy scale of new physics. We assnmed 
the same (gange) symmetries for the as for Cq. This implies that the d^ term (and 
all terms with an odd d) in (j^) are absent. If we fnrther assnme that the gj are of the same 
order of magnitnde as the standard model conplings g,g' and e and compare the Lagrangian 
with the one defining the a-parameters |^, we read off the order of magnitnde for the 
a-parameters, 

aw:Oiw^,(y.B^ = O [m‘^/ hdj . (7) 

Assnming A > 2 TeV (and conseqnently restricting onrselves to scattering energies np to 
Mv^Va < 2 TeV), the order of magnitnde for the a-parameters is 




< 


o (io-3). 


( 8 ) 


The restrictions derived from partial wave nnitarity applied to vector boson scattering 


amplitndes are 


saw 

Ml 


< 



saw^ 

~W 


S 15.5, 


sa 

'W 


< 49 


(9) 


where we have introdnced s = My^y^. For a/s < 2 TeV the nnitarity bonnds (|) are 


|o^w| Z 0.031, Z 0.025, |q^b<i>| Z 0.079. 


( 10 ) 


These limits are larger than the valnes in Eq. (H) for the a's which we expect from the 
effective Lagrangian ansatz. Therefore, if the conplings are not larger than expected from 
the effective Lagrangian ansatz, nnitarity is not violated for energies a/s < 2 TeV. In 
O, 1^, ^ 1^ a form factor assnmption is made in order to avoid violation of nnitarity. 


In onr fits we follow the simple prescription to vary the conpling parameters within their 
nnitarity limits only. In fact it will tnrn ont that within the 95% CL limits the nnitarity 
limits are never reached. Tims, in order to derive sensible experimental bonnds on the 
anomalons conplings, one does not have to nse form factors for which additional (nnknown) 
parameters mnst be introdnced. 

If one nevertheless introdnces a form factor, the conplings ai which are to be inserted 
in the expressions for the cross sections are energy dependent. They are related to bare 
(energy independent) conpling constants, by 


a,- 


di = 


1 + ^ 


FF 


n ■ 


( 11 ) 


The bare conpling constants are those which appear in the Lagrangian. A nsnal choice for 
the exponent n in (pTf ) is n = 2. Similar to A in (^ Kpp is an energy scale for new physics. 
The nnitarity limits for the parameters are obtained by inserting (pT]) into (^. We nse 




































n = 2 and minimize the maximum value for |q!°| with respect to s. The minimum occurs at 
s = h?pp and the unitarity limits are given by 

|a^| < 7Q{M^/Alp) ~ 0.123, < 62(M^/A|^) ~ 0.100, 

< 196{M^/Alp) ~ 0.316. (12) 

The numerical values in (^) are for App = 2 TeV. At multi-TeV colliders the cross section 
for hxed 7 ^ 0 is very different from the cross section for hxed a* and the obtainable 
bounds on the a* are very much tighter than those for the a°. The distinction between the 
two models does however not very much affect the analysis of present Tevatron data since 
there the form factor is close to the value 1 as a/s can hardly be greater than 0.5 TeV. 

2.4 Results with Anomalous Couplings 

Fig. 0 shows the comparison of ^/-annihilation and vector boson fusion in the presence of 
anomalous couplings. We show the results for the relevant processes of W^Z and 
production and for W^W~ production. In addition, we present a plot for pro¬ 

duction. We sum over the charge conjugated hnal states i.e. discuss the cross sections for 
W^V = W+V + W-V and = W+W+ +WW- production. We have also summed 

over all V 1 V 2 pairs. We only vary one coupling at a time. Only those couplings which lead 
to enhanced terms at high energies (i.e. of O(ais) or O(afs^)) in the qq' cross section are 
varied. Varying the other couplings leaves the qq' cross sections virtually unchanged. For 
production we vary all couplings. We choose a single non-zero magnitude for each 
of the couplings which is already quite large for the effective Lagrangian expectation, (|), 
but which is still below the unitarity limit (|^). For aw and aw^ we take \ai\ = 0.01. For 
aB<s> we take |as$| = 0.03. For the relevant processes of hF^V production, we choose a 
negative and a positive value for the coupling if there is an enhanced term linear in the 
coupling. 

The main conclusion from Fig. ^ is that vector boson scattering is only marginally 
important even if the anomalous couplings are different from zero. When constraining 
anomalous couplings using these processes, vector boson scattering might therefore well 
be omitted. The non-enhanced terms (a^^ in IV^Z-production, aw^ and ap^ in PF^y- 
production) are unlikely to lead to any observable effect at the LHC. Fig. (d) shows that 
the effect of anomalous couplings for like-charge PF^-pair production is not very large. 


3 Parameter Fits for Anomalous Couplings 


In this section we present parameter hts to hctitious standard model data and derive limits 
for the anomalous couplings. Refering to the conclusion of Section Q, we will take into 
account only the contribution from qq' annihilation. First we consider PF^y and PF^Z 
production separately. These are the experimentally relevant production processes |]T3 . 
The detection of a PF+PF“ pair is experimentally plagued by a large background of tt 
production with the subsequent decay of a top quark into a PF^ boson and a b quark [jlS . 
We use the general parametrization of the triple gauge boson vertices [^, ^ in terms 

of seven free parameters, thus allowing for C- and P-violation. Then we present a £t to 
combined PF^y and PF^Z “data” for the three parameter gauge invariant model. We take 
into account the full correlations among the parameters. Before we proceed we present the 
unitarity limits for the set of couplings which we are using |^, ^ . As far as we know, these 
limits have never been given before. 

















PP ^ ^ 3^4 + X at Vs = 14 TeV 



0 0.5 1 1.5 2 0 0.5 1 1.5 2 


Mv3vjTeV] 

Figure 7: The cross-sections for W^Z{= IF+Z + W-Z),W^-f,W+W- and W^W^{= 
W~^W~^ + W~W~) prodnction as functions of the invariant mass My^Vi of fhe prodnced 
vector boson pair for pp-collisions at y/s = 14 TeV. Various values of the anomalous couplings 
have been chosen. Separately shown are the contribntion from gg'-annihilation and the 
(snmmed) contribntion from vector boson fusion processes. A psendorapidity cnt of ?7 = 1.5 
has been applied. 
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3.1 Unitarity Limits for dy, xy,yy, zy, zi 
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Theoretical bounds on anomalous couplings can be obtained by applying partial wave uni¬ 
tarity to the amplitudes for qq' V 3 V 4 . Inequalities derived from the requirement of partial 
wave unitarity have been given in |^. The inequalities have been written in terms of “re¬ 
duced amplitudes” for qq' —>• and qq' scattering. The reduced amplitudes 

have been given in terms of the parameters gY, i^v, dX, 9^, and Ay, where V = j 
or Z. By comparison of the Lagrangians of |^| and we hnd the following equivalence 


between this set of parameters and the one we are using: 


9^-1 = 


6 


V 


nSM’ 

9v 

J 


9i 


Ky = 


Z 


IV 


nSM^ 

9v 

J 


z. 


2V 


9^ 


Ay — 


9^ = 


Vv 

nSM^ 

9v 

Zv 


KV — 1 — 


6y -|- Xy 


MY 


^3V 


(F 


9Y^MY. 

^ 


+ i 




(p2 


„SAr 

9v 




„SM 

9v 


) . (F^ 




9^ 


MY. 


— i 


M^y*) 


9r 


MYv 


Ay = 


-2 


9v 


(13) 


In (0), F\ My/* and My are the squared invariant masses of the IT^, the W and the V, 
respectively, entering the trilinear vertex and = 1, = {cw/sw)- Because the two 

parameter sets are only equivalent up to possible form factors^, the equations (^) contain 
the kinematic variables Myr* and We checked that with the replacements (|^ the 
expressions for the amplitudes for W^Z and W^W~ production in terms of the two 

sets translate in the correct way. 

The following table summarizes the symmetry properties of the parameters under C and 
P transformations: 


^V,Xy,yy 

Zy 

z'lv 

^3V 

C,P 

(CP) 

iZ,P 

c,p 


If electromagnetic gauge invariance is demanded the following parameters vanish, 


I'll). 


5.y = 0 , 


2 )^ = 0 . 


(14) 


Assuming that only one anomalous coupling at a time is different from zero we extract 
the unitarity bounds shown in Tables |l| and ^ from the bounds on the reduced amplitudes 
For W^Z production we neglected terms of 0{MY//s). For the form factor case we 


m 


used (pd]) with n = 2 and minimized the unitarity bounds with respect to s. However, for 
z!^, zY and {z'^z)^ ^^e value of s at the minimum is greater than A^p. For these cases we 
quote the unitarity limit for s = A^p- The bounds shown in Tables IT] and 0 are weaker than 


those derived from vector boson scattering because in the latter processes the amplitude is 
in general quadratic in the couplings while for qq' —>■ V3V4 it is at most linear. 


3.2 Present Direct Limits 


At present, direct limits on the couplings have been obtained by the CDF and DO collabo¬ 
rations at Tevatron |Q, Q and by the LEP 2 collaborations [^, Q . Table ^ summarizes 
the most stringent bounds which were attained at the Tevatron. 


^Form factors can be introduced by adding terms with two or a larger even number of derivatives on the 
fields to a Lagrangian with constant couplings. These terms are equal to a power of a squared invariant mass 
(or even a product of powers of several squared invariant masses) times the interaction term of the original 
Lagrangian. In order to compare the interaction terms of the Lagrangians of and ||^ the terms of one 
Lagrangian have to be re-grouped (by using partial integrations and tensor identities). Two derivatives on 
a field appear in some of the re-grouped terms. This introduces the P^, and My dependences in 



































Para¬ 

meter 

Unitarity 

limit 

v^ = 

2 TeV 

Para¬ 

meter 

Unitarity 

limit 

KpF — 

2 TeV 

\x^\ 

I1/7I 

\z^\ 

\4^\ 

K/{y/2p) ~ 6.0/v/i, S > 

V2K/p ~ 12.0/v/s 

y/2K Mw / — 0.96/s 

y2ir/[z/((s/M^) - 1)] ~ 0.077/s3/2 

\f2Klv ~ 12.0/v/s 

v^77/z/ ~ 12.0/v/i 

3.0 

6.0 

0.24 

0.0096 

6.0 

6.0 

UOI 

U7I 

ll/°l 

U7I 

l(4)°l 

I( 47 )°I 

96/(3y3Ai.j.) 

192/(3y3Airi.) 

3.84/A^^ 

0.308/A^^ 

192/(3y3Airi.) 

192/(3y3Airi.) 

9.2 

18.5 

0.96 

0.039 

18.5 

18.5 


Table 1: Unitarity limits for '-yW^W couplings without and with a form factor derived 
from partial wave unitarity for qq' —IT^y. K = v^sv^^/(q; ^l-Ml./s), u = ^s/M^ + 1. 


Para¬ 

meter 

Unitarity 

limit 

v^ = 

2 TeV 

Para¬ 

meter 

Unitarity 

limit 

ApF — 

2 TeV 

\Sz\ 

{2K/sw){M^Is) ~ 1.54/s 

0.39 


6.16/A|^ 

1.54 

i^zj 

^/2K{cw / sw){Mw / \/s) — 12.0/ys 

6.0 

uo 1 

192/(343Ai.j.) 

18.5 

\yz\ 

^/2K{cw/sw){M'^/s) ~ 0.96/s 

0.24 

bll 

3.84/A|^ 

0.96 

jzzj 

V2K{cw/sw){M^/s^/^) ~ 0.077/s3/2 

0.0096 

141 

0.308/A|^ 

0.039 

kiz 

{2K/sw){M^/s)-IM/s 

0.39 

l(^4)°l 

6.16/A|,^ 

1.54 

\^2z\ 

y/2K{cw/sw){Mw/y/s) ~ 12.0/ys 

6.0 

l(44°l 

192/(343Ai.^) 

18.5 

ksz 

K/{V2sw){M^/s^/^) ~ 0.044/s3/2 

0.0055 

l(44°l 

0.176/A|^ 

0.022 


Table 2: Unitarity limits for ZW^W couplings without and with a form factor derived 
from partial wave unitarity for qq' W^Z. K = /{acw)- 
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diboson pair 

W-f 

W+W-/WZ W+W- 

W+W-/WZ 

assumptions 

none 

Sz = 0,xz = l^ 


y-y = yz = o,x^ = 0 

results 

-1.4 < < 1.4 

—0.5 < < 0.5 

—0.4 < Xj < 0.6 
-0.4 <y^ < 0.4 

—0.9 < x^ < 1.0 
-0.7 <y^< 0.7 

- 2.5 <6z < 2.7 
-4 < Xz < 4 

ApF/TeV 

1.5 

2 

2 

1 


Table 3: Results from [Q for the 95% CL limits on anomalous couplings obtained from 
two-parameter fits to data of diboson production in pp collisions at y/s = 1.8 TeV. The 
bounds take into account possible correlations between the two htted parameters. The C- 
or P-violating couplings were assumed to be zero. The value of which was used in the 
hts is indicated in the bottom row. 


The LEP 2 collaborations recently gave a preliminary limit for 

-0.3 < aw# < 0.4, 95% CL, 


where a^^ = aw = 0 was assumed. Adopting a two-parameter model [^] which is equiva¬ 
lent to aw#, as$ % 0 and aw = 0, the following preliminary limits were obtained at LEP 2 


5e 


\5z\ < 1.9, 


95% CL. 


—2.5 < < 3.8, 

These limits take into account the correlations between the two parameters. No form factor 
was used. 

For the three parameter 


The hnal sensitivity of LEP 2 has been estimated in 37, 


gauge invariant model the following result was obtained for a run at ^/s = 190 GeV with 
an integrated luminosity oi C = 500 pb“^ [p7|: 


—0.20 < aw < 0.24, —0.19 < aw# < 0.13, —0.35 < as^ < 1.05. 

These bounds are at 95% CL and take into account all correlations. 


3.3 Fitting Procedure 

We performed hts to the Mw and pr distributions of the cross sections, where pr = q sin 6 
is the transverse momentum of a produced vector boson. The pt distribution was calculated 
according to 


da 

dpT 


(hiha ^ qq' V3V4, Sfe/j)|cut 


Pt J zt 


2 min[-(l/2) lna;,yobT)] 




X 


WLa.x[{l/ 2 )\nx-yo{zT)] 

{qq' hsiA, zt) + - — n{qq' hsiA, —zt) 


dcosO 


dcosO 


h2 

(15) 


In ([IsD , Zt = Y^l — {pT/qY is the magnitude of cos 6 for the given pt and Xmin and Xmax are 
determined by 


Xmia = max 


(Ms -F M^f/shh, x{pt] 



















































pp ^ W=^V + X at Vs = 14 TeV 



Figure 8: The cross sections for pp —>■ 1F^7 + X (a) and pp —>■ W^Z + X (b) as a function 
of the transverse momentum px of a produced vector boson in the standard model and for 
various values of the anomalous couplings. The cross sections have been multiplied by the 
branching ratios for the decays of massive vector bosons into two generations of leptons. 


^max 


min 


(2 TeVjVskk 


( 16 ) 


In (^) we included the upper bound of 2 TeV for -y/s. Vq^zt) in ([l5|) is determined by the 
pseudorapidity cut. In the high-energy limit it is given by 

Vo^Zt) ~ 7 - tanh“^(2:r). (17) 


The function x{q^) in ([I^) is given by 


x{q^) = — - ^^(18) 

Shh 


and -dmin in ([IGD is determined by the relation tanh(p) = cos'dmin- Q is the variable defined 


m 


We assumed that the hF^, Z particles are identified by their decays into two generations 
of leptons each. We used the following branching ratios, 


10 . 


Z^l+l- 3.37% 


We used no other cut than p = 1.5 on the produced bosons. Figs. || (a) and (b) show the 
cross sections for pp —>• + X and pp —> W^Z + X, respectively, multiplied by the 

branching ratios as a function of pt in the standard model and for various values of the 
anomalous couplings. No form factor was used. 
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Bin Nr. 

1 

2 

3 

4 

5 

6 

7 

Mw [TeV] 

1-2 

0.8-1 

0.7-0.8 

0.6-0.7 

0.5-0.6 

0.4-0.5 

0.3-0.4 

Nsm 

C ■ Br ■ a 

W±7 

95 

126 

134 

247 

499 

1154 

3320 

W^Z 

21 

29 

31 

59 

121 

285 

823 

zz 

3.3 

4.7 

5.2 

9.8 

20.4 

49 

147 

Z 7 

33 

45 

48 

89 

182 

426 

1247 

77 

214 

287 

304 

556 

1110 

2518 

6965 


Bin Nr. 

1 

2 

3 

4 

5 

6 

7 

8 

VT [TeV] 

0.5-1 

0.4-0.5 

0.35-0.4 

0.3-0.35 

0.25-0.3 

0.2-0.25 

0.15-0.2 

0.1-0.15 

Nsm 

VP±7 

32 

51 

55 

103 

208 

473 

1275 

4603 

W^Z 

7.6 

12 

13 

24 

48 

108 

275 

821 

ZZ 

1.6 

2.5 

2.7 

4.9 

9.8 

21.5 

53 

151 

Zl 

17 

25 

27 

50 

101 

225 

595 

2048 

77 

116 

172 

187 

345 

698 

1600 

4460 

17600 


Table 4: The numbers of standard events in 7 bins over the invariant mass Mw and in 8 bins 
over the transverse momentum pt for the processes pp —> V3V4 + X sX y/s = 14 TeV with 
a cut \p\ < 1.5 on the pseudorapidity of the produced vector bosons and the requirement 
^/s < 2 TeV for the pt distributions. An integrated luminosity of £ = 10^ pb~^ has been 
assumed. For massive vector bosons a decay into two generations of leptons was assumed. 
All results were obtained in the Born approximation. 


To estimate the number of events at the LHC we assume an integrated luminosity of £ = 
10^ pb“^. We arrange hctitious standard model data into bins. For the Myv distribution 
for pp —> + X we hnd that there is less than 1 event for Myv > 2 TeV and ~ 21 events 

in the interval 1 TeV < Myy < 2 TeV. We choose this interval to be the hrst bin. The 
other bins and numbers of SM events for and production are shown in Table 

1^. We also show the numbers of events for ZZ, Z7 and 77 production^. The accuracy of 
the numbers due to numerical integration is 1%. We proceed to arrange the data for the pt 
distributions into bins. Since pt — Myy/2 for scattering at right angles and large invariant 
masses we choose the limits for the px bins equal to half the limits of the Mw bins. In 
addition we dehne an eighth bin. The results are also shown in Table 

To calculate the non-standard effects we wrote the number of events in each bin as a 
power series in the anomalous couplings, 

N{as) = Nsm + atNi + Y ( 19 ) 

i i,j 

where the a* are the anomalous couplings. Table ^ shows the coefficients for the C- and 
F-conserving couplings and for z.y in a bin of the px distribution comprising 0.4 TeV < px < 
1 TeV. 

In each bin we calculate 

Ax^ = -2\n(^-^y ( 20 ) 

where C is the likelihood function for the data in this bin, assuming that a theory with 
particular (non-zero) values for the anomalous couplings is the correct theory, and Cq is 


®We neglected potential contributions from gluon fusion 




































IT±q 

SM 

£7 

Vz 

^7 

q2 

x^ 

y'l 

^7 

X'fjj'Y 

X.yZ.y 

jj'yZ^ 

82 

-2.5 

0 

-6.1 

63.5 

233 

497 

31.8 

0.1 

0.01 

IT±Z 

SM 

Sz 

Xz 

yz 

A2 

^Z 

y2 

Xz 

fz 

SzXz 

^zVz 

xzVz 

20 

-12.7 

-3.0 

-0.68 

6.3 

4.1 

15.3 

2.3 

0.57 

2.0 


Table 5: The coefficients iVj and Nij, defined in Eq. (|T^), for the numbers of produced 
and W^Z pairs in a bin with 0.4 TeV < < 1 TeV for the rescaled coupling parameters 

■ 10, ■ 100, ■ 10^ and Sz = Sz ■ 100, xz = xz ■ 10, yz = Hz ■ 100. The 

numbers of standard events are also shown. Sample usage: For y^ = 2 ■ 10“^ y^ = 0.2 
and all other anomalous couplings equal to zero there are 82 + 233 ■ (0.2)^ ~ 91 events 
in the bin, corresponding to one standard deviation from the standard model. 

the same function assuming that the standard model is correct. is a measure of the 

probability that this particular model can still describe the (standard) data. If the number 
of events (in the bin) is greater than 50, we calculate C according to a Poisson distribution 
of the total number of events. 


< N 

C = PN = ^^e-<^>. (21) 

In (0), < N >= iV(a's) is the number of events predicted by the particular non-standard 
theory and N = Nsm is the number of standard events (=the number of “measured” 
events). If the number of events is smaller than 50 we generate the N events in the bin, 
i.e. we calculate the phase space points fli, i = 1... N, a.t which the standard events would 
be located in the bin. hi represents Myv or pt for the two distributions, respectively. We 
then use the method of extended maximum likelihood (EML) to calculate C. The likelihood 
function of the EML is given by 


N 


^EML — PN np(s^..'3). 


( 22 ) 


with ptv from (^Tj) and p{^li, a) is the probability of finding the ith event at the phase space 
point hlj, assuming that the theory with the parameters a is correct. p{Qi,d) is given in 
terms of the differential cross section by 


p(fli,a) =-^(fli), with cr=/' 
(J I Jhi 


da 


bin cir2 


dhl. 


(23) 


where a and da/dVL are evaluated in the non-standard theory. 


3.4 Results 

Fig. ^ shows the projections of the = 1 and Ay^ = 4 confidence regions on the 

x..), = 0, ?/..y = 0 and z.y = Q parameter planes as the result of a three parameter fit to the pt 
distribution of pp hP^q -|- X. The parameters ^27 and z'.^^ were set equal to zero. As 
the parameters are uncorrelated, the projections are equal to the sections of the confidence 
regions with the planes. If the Myy distribution is used instead, the regions expand by 
a factor of 1.1 to 1.15 in each dimension. If only the total number of events in each bin 
is subjected to the fit (instead of using the EML method), the regions expand by a factor 





















of about 1.2 in each dimension. If a four parameter fit to x^,y^,z^ and Z 2 ^ is performed 
instead, the projections and sections stay the same as in Fig. y. The parameter z'^^ does 
not contribute to production. We assumed 6^ = z[^ = Q because of electromagnetic 

gauge invariance. 

Fig.0 shows the projections of the conhdence regions on the parameter planes 5z = 0, 
xz = ^ and yz = ^ and the sections of the regions with the planes as the result of a three 
parameter £t to the px distribution of pp —>• W^Z + X. The parameters and 

were set equal to zero. The figure displays the correlations among the parameters. As a 
result of the correlations the sections are smaller than the projections. A four parameter £t 
which includes yields identical results. If a seven parameter £t (including also 
and is performed, the projected conhdence region expands by no more than 4% in any 
direction except for the positive 5^ direction for = 1 where it expands by ~ 30%. 

Figure ^ shows the projections and sections on the aw = 0, aw# = 0 and a^^ = 0 
planes from a simultaneous £t of the px distributions of pp +X and pp PF^y + X 

to the three parameter gauge invariant model. The unitarity limits for the parameters are 
also shown. The conhdence regions he inside the unitarity limits. 

The use of diherent parton distribution functions leads to small theoretical uncertainties 
(< 1%) in the conhdence regions. These uncertainties could be reduced by subjecting ratios 
of cross sections, e.g. (j(PF^l/)/(j(yy), to the ht. Due to the additional statistical error 
induced by the reference cross section (i.e. o'(yy) in our example) the conhdence levels 
derived from the ratios are, however, several tens of percent wider than those derived from 
the absolute values of the cross sections. We do not use ratios for our hts. 

We repeat our analyses using a form factor. We project the conhdence regions on the 
parameter axes. This results in 95% (for = 4) and 68% (y^ = 1) conhdence limits for 
the parameters. Table ^ summarizes our results. If we repeat the hts for y = 3 the bounds 
are only slightly ahected: the diherences between the maximal and minimal values of the 
couplings change by at most 20% compared to Table ^ In general the diherences decrease. 

The 95% conhdence limits which we obtain for the alternative set of parameters A^ff, Xkz 
and of (instead of Sz,xz and yz) are: 


0.0028 < Ac/f < 0.0080, 
-0.062 < Akz < 0.044, 
-0.0041 < Az < 0.0041, 


-0.0052 < A^f’° < 0.024, 
-0.13 < A4 < 0.062, 
-0.0103 < A| < 0.0104. 


(24) 


We compare our results with previous investigations. Sensitivity limits achievable at 

Fits to the px distribution of 




the LHC were previously presented in 
hctitious data for pp —> W^Z + X at a/s = 14 TeV were performed in [^. The 95% CL 


limits presented there, using a form factor with Kpp = 3 TeV and n = 2 and based on the 
Born level prediction were 

-0.0048 < A^f’° < 0.0164, -0.120 < Ak^ < 0.092, -0.0082 < A^ < 0.0084. 

If we repeat our three-parameter ht with App = ^ TeV we obtain a similar result, namelyQ 
-0.0039 < A^f’° < 0.0140, -0.090 < A4 < 0.053, -0.0067 < A| < 0.0068. 


An SSC analysis using a form factor for pp 
repeat our analysis with the parameters used in 


PF+y -P X can be found in |2^. If we 
(y^ = 40 TeV, PF’*' decays to only one 


^ Different cuts were used in |^. In particular, a pseudorapidity cut of 77 = 3 was applied on the decay 
products of the vector bosons. If we repeat our analysis for 77 = 3 and include, as in only W^Z 
production our limits change by less than 3%. 






















Ax^ —surface, projections 
from pp ^ W 7 , 7] = 1.5 


= 0 plane ^ 0 plane 



= 0 plane 

* 10 ^ 



Figure 9: The projections of the = 4 and Ax^ = 1 conhdence regions on the and 
7 parameter planes from a three parameter £t of the pt distribution of pp tF^7 + X 
at \/s = 14 TeV with a cut of 77 = 1.5. An integrated luminosity of A = 10® pb“^ and a 
leptonic decay of the tF^ boson into two generations of fermions was assumed. All other 
parameters were assumed to be equal to zero. 


21 




Ax^ —surface, projections and sections 
from pp ^ WZ, 7] = 1.5 


^ 2^0 plane 


cD^ ^ 0 plane 




Figure 10: The projections of the Ax^ = 4 and Ax^ = 1 conhdence regions on the Sz,xz 
and yz parameter planes and the sections of the regions with these planes from a three 
parameter £t of the pt distribution of pp + X at ^/s = 14 TeV with a cut of 

p = 1.5. An integrated luminosity of A = 10^ pb~^ and a leptonic decay of the tF^,Z 
bosons into two generations of fermions was assumed. All other parameters were assumed 
to be equal to zero. The sections are drawn in the same way as the projections. They can 
be distinguished from the projections as they always he inside them. 
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—surface, projections and sections 
from pp ^ WZ and pp ^ W 7 , 77 = 1.5 

= 0 plane = 0 plane 

*100 *100 



= 0 plane 

*10 



Figure 11: The projections of the Ax^ = 4 and Ax^ = 1 conhdence regions on the avn, o:w^ 
and as^ parameter planes and the sections of the regions with these planes from a simul¬ 
taneous £t of the pt distributions of pp -|- X and pp tF^7 -|- A to the three 

parameter gauge invariant model for ^/s = 14 TeV and a cut of 77 = 1.5. The unitarity 
limits are also shown. An integrated luminosity of £ = 10^ pb~^ and a leptonic decay of 
the tF^, Z bosons into two generations of fermions was assumed. All other parameters were 
assumed to be equal to zero. The sections are drawn in the same way as the projections. 
They can be distinguished from the projections as they always he inside them. 
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no form factor 



with form factor 


parameter 

68% CL 

95% CL 

parameter 

68% CL 

95% CL 

min 

max 

min 

max 

min 

max 

min 

max 



5z ■ 100 

-0.24 

0.87 

-0.51 

1.48 

■ 100 

-0.49 

3.40 

-0.99 

4.41 

Xz ■ 10 

-0.68 

0.38 

-1.26 

0.82 

x% ■ 10 

-1.77 

0.64 

-2.76 

1.29 

yz ■ 100 

-0.46 

0.46 

-0.75 

0.76 

I/O ■ 100 

-1.44 

1.43 

-1.93 

1.94 

;^z ■ 10=^ 

-0.26 

0.26 

-0.45 

0.45 

z| ■ lOO 

-0.85 

0.84 

-1.30 

1.29 

■ 100 

-0.74 

0.74 

-1.21 

1.21 

(zjz)° ■ 100 

-2.8 

2.8 

-3.1 

3.1 

^2Z ■ 10 

-1.35 

1.35 

-1.79 

1.79 

(^2z)° ■ 10 

-3.0 

3.0 

-3.4 

3.4 

4z ■ 10" 

-1.47 

1.47 

-2.55 

2.55 

(4z)° ■ 10" 

-4.9 

4.9 

-7.4 

7.4 

x^ ■ 10 

-0.23 



■ig 

a;0-10 

Him 

0.56 

HH 

0.77 

y^ ■ 100 

-0.131 


■B 

wm 

|/° ■ 100 

BB 

0.34 

BB 

0.50 

■ 10^ 

-0.075 

OBI 

-0.129 

0.135 

zO ■ lOO 

BB 

0.25 

-0.37 

0.39 

4-v ■ 10 

-0.30 

BB 

-0.44 

0.44 

(z',)0 ■ 10 

-0.45 

0.45 

-0.63 

0.63 

aw ■ 100 

BSSi 

nnra 

BSB 

0.20 

a^ ■ 100 


0.31 

bb 

0.47 

aw^ ■ 100 



BB 

0.26 

a'^^ ■ 100 

^BB 

0.22 

BB 

0.42 

q;b<J) ■ 10 

-0.22 

BB 

Bn 

0.48 

■ 10 

bob 

0.53 

-0.47 

0.74 


Table 6; The projections of the Ax^ = 1 (68% CL) and Ax^ = 4 (95% CL) conhdence 
regions on the parameter axes as the resnlts of a seven parameter £t of the pt distribntion of 
pp hh^Z + X to Sz, xz, Vz, zz, x'^z^ ^2Z ^ parameter fit of the pt distribntion 

of pp hL^7 + X to and ^2^ and a three parameter £t of the combined pr 

distribntions of pp W^Z + X and pp fT^7 + X to aw, aw<s> and aB<s>- The form factor 
results are for AiT’iT’ = 2 TeV and n = 2. 
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lepton family, C = 10^ pb“^ and using only the information from the total number of events 
in each bin) and use a cut of r/ = 2.5 we obtain -0.17 < < 0.21 and -0.021 < A° < 0.020 

at 95% CL0. These bounds are tighter by a factor of 1.5 to 2 than the ones obtained in p3 


In 


an LHC bound on was derived, assuming = ^. This bound was derived 


from the 0{as) prediction for the cross section, but from Table IV of |24 


we deduce that 
< 0.06PI ("we do 


the la bound which would be obtained from the Born approximation is \x^ 
not use a form factor for this comparison). This bound is wider by a factor of in between 


two and three than our bound. This can be explained by the fact that in the assumed 
luminosity was only C = 3 ■ 10"^ pb~\ only IT+y production was considered, the fitting 
procedure was simpler (only one bin was taken) and different cuts were used. Including the 
0{as) corrections reduces the sensitivity to x^ by about a factor of two 

The limits which were derived in are much larger than ours because the discovery 
criterion employed there is much stronger than ours. We note that the chiral Lagrangian 
parameters Xg and Xg used in 


10 


are 


identical to and 05$, respectively!^ The explicit 


connection is given by 


a 


a 


Stts 


—Xg — —aw^i 


w 


Stts^ 


Xg = -Q!S$. 


(25) 


Conclusion 

We showed that the rate for vector boson fusion production of vector boson pairs at the 
LHC is at the order of 10% to 20% of quark antiquark annihilation production and might 
thus be neglected in an estimate of the pair production cross sections. This result was 
obtained by applying an improved formulation of the effective vector boson approximation 
(EVBA). It agrees with the result of a calculation in which the complete set of diagrams 
was evaluated instead of performing an EVBA. Previous calculations in which the EVBA 
in leading logarithmic approximation was used overestimate the contribution from vector 
boson fusion by a factor of 3. 

We derived conhdence intervals for the full set of anomalous IV^IV“7 and W~^W~Z 
couplings, including C- and P-violating couplings, from fits to the standard IV^y and W^Z 
production rates expected for the LHC. In addition we derived conhdence intervals for a 
three parameter SU{2)l x U{1)y gauge invariant dimension-six extension of the standard 
model. We performed multi-parameter hts in which the full number of anomalous couplings 
was varied. Our limits thus take into account all possible correlations among the effects of 
the various possible couplings. We derived limits with and without making a form factor 
assumption. We compare the limits with the unitarity limits for the production of vector 
boson pairs with invariant masses smaller than ^/s = 2 TeV. It turns out that all 95% 
conhdence limits lie inside the unitarity limits whether a form factor is used or not. It is 
therefore not necessary to use a form factor in order to avoid violation of unitarity. The 
limits which we obtain without using a form factor are a factor of 10 (for x^ or as^) and 100 
(for {6z,y-y) or {aw<s>,o:w)) stronger than the present experimental limits or limits which 
can be attained at LEP 2. 

Adopting an ehective Lagrangian approach, together with an assumption about the 
energy scale at which new physics occurs, provides us with an order of magnitude estimate 
for the parameters a-w, aw<i> and We hnd that for aw and aw^, the limits which can 
be obtained at the LHC are of the same order of magnitude as this estimate. It might thus 

®We required pr > 200 GeV. 

® We assumed a quadratic dependence of the number of predicted non-standard events on this coupling. 

^°This is true as far as only the trilinear vector couplings are concerned. 




























be possible to observe non-zero values of these coupling parameters, should they exist, at 
the LHC. 

In appendices we give analytical expressions for cross sections for vector boson pair 
production with anomalous couplings for qq' annihilation and vector boson fusion. Our 
expressions manifestly show the effects of the couplings at large scattering energies. 
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A Cross-Sections for gg'-Annihilation 


The standard model differential cross sections to O(a^) for qq' —>■ qq' —>■ PF^y, 

qq PF^1F“ and qq ^ ZZ have been hrst given in [^[5 ^ form in which good 

high energy behavior is manifest and including the aiy-interaction, all cross sections for 
qq' —> V 3 V 4 can be found in [^]. For arbitrary vector boson self-interactions all cross 
sections and helicity amplitudes have been recently given in [H 


We give here the formulas for the differential cross sections for the qq' processes which 
receive contributions from anomalous vector boson self-interactions, qq' —>■ PF^Z, qq' —>■ 
PF ^7 and qq PF^PF“, in a form in which the high energy behavior is manifest. As in 
[^, we have explicitly carried out the high energy cancellations among different diagrams, 
also (as far as possible) for the non-standard terms. We use the general C- and P-conserving 
vector boson self-interactions compatible with Lorentz-invariance and electromagnetic gauge 
invariance in terms of the parameters x^,y^, 6 z,xz and yz- In addition, we include the 
contributions from and 2:270 for PF^y production and the contributions from zz, z'^z^ 
and z'j^z fo^ PF^Z production. The differential cross sections for qq ZZ, yZ, yy can be 


found in 30 


For q^ —>• PF^F, V = y, Z, the cross sections contain an overall factor of where 

Vqqi is the element of the CKM matrix for the mixing of the quarks q and q'. For qq —> 
PF’''PF“, the quarks have to be of the same flavor. We give the cross sections averaged over 
colors and spins of the initial quarks and summed over the helicities of the final state vector 
bosons. The cross sections given here agree with the expressions given in and 


We denote the left- and right-handed couplings of the Z-boson to the quarks by 


We also use the symbols 


Lu — ^ 3 '^'^’ 

2 2 

Ld — —^ + 


x" — 1 — —1 ■ 

'3 ~ '3 ~ -‘-5 




W 


Rd ^ xsi 


2 1 
Qu g) Qd g- 


(26) 

(27) 

(28) 


^^Also the a;.y-terins for the TF^y production cross section have been given there. 
^^The parameter does not contribute. 

After correction of misprints. 


















The Mandelstam variable t will be defined below for each process and u is defined by 
s + t + u = M| + M|. The scattering angle 9 is the angle between the three-momenta of 
the two particles which define t. We further use the variables 


1 


2 (M| + M|) (M| - Mi )2 

I o • 


and rj, where rj = for production. We treat the processes qq' 

qq' —>■ together because similar functions are involved. 

The differential cross sections are given by: 


and 


qq 


W^Vi 


1. qq' —» 


_ ^^ 11 / 12/92 f 1 
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tu 


w 2 
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SSw 
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(0) + ‘^x^y^Axy{9i) -f- -f (^ 27 )^) ^* 2 ( 0 ) 


>, (29) 


where we have defined t as 


t = 


(Pu-Pw+Y 
{Pu-Pw-Y ’ 


(30) 


for the two charge conjugated processes, respectively. In, (|30|) Pi denotes the four-momentum 
of the particle i. 

2. qq' 
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t = [ . (32) 

\ {Pu-Pw-) 

The invariant functions for qq' —> f+^V 4 for the terms linear in the anomalous couplings 
are given by 
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and the functions for the terms quadratic in the couplings are given by 
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qq W+W- 
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with 


The invariant functions for qq 
are given by 
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(36) 
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and the functions for the terms quadratic in the couplings are given by 
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Leading Terms in qq' —>■ W^Z 

Leading Terms in AI p^, qq' —>■ 
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Table 7: The leading behavior for s S> of the helicity summed squared amplitude 
|Afp for qq' —> W^Z and qq' —> Shown is the leading behavior of the different 

terms in proportional to the different combinations of the anomalous couplings. If 

a term has a different sign for W^{ri = —1) and W~{ri = 1) production this is indicated 
by the factor r]. The terms can be even or odd in cos 6. If they are odd this is indicated 
by the factor cos6'. A slash indicates that a term is not present. Sample usage: \Ai\^ = 
0 ( 5 °) + dzOjs/Ml^) + b\0[s ^!for q^ —» W^Z if only bz is non-zero. 


-6x2 

Bxy 

By2 



(38) 


Table ^ shows the behavior for s S> for the helicity amplitudes for qq' —> W^V. We 
note that the terms which are proportional to cos 9 give no contribution to either the pt or 
the Mvv distributions. 


B Cross-Sections for W^V2 —> Scattering 


We give expressions for cross-sections for V 1 V 2 —> V 3 V 4 in a high-energy approximation (to be 
described below). We restrict ourselves to the relevant processes oiWZ and IT 7 production 
(in the following we simply write W instead of W^). Thus, we only give the cross-sections 
forlTZ —> ITZ, IT 7 —> WZ^WZ IT 7 and IT 7 —> IT 7 . Helicity amplitudes for processes 
V 3 V 4 in the high-energy limit of the GIDS model can be found in [^, ^ ^ 




They 


have been obtained from the exact Born-level amplitudes by an asymptotic expansion for 
s 3> where s is the scattering energy squared. The expansion has been carried out at 
a fixed scattering angle 6 . Therefore, also \M‘^/t\ and \M‘^/u\ have to be small parameters. 
We note that for scattering energies s > 0.8 TeV the parameters \M‘^/t\ and \M‘^/u\ are 
smaller than 0.2 for all scattering angles if a pseudorapidity cut of 7 < 1.5 is applied. 

Since we assume that the couplings are small, at = 0(M^/A^), we only keep those 
anomalous terms in which each power of an anomalous coupling is enhanced by a factor of 
For this purpose we define the parameters 


s/4M2 
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dj — 


4M^ 




(39) 


The assumption is equivalent to assuming a* = 0(1) or smaller (since 

s < A^). In addition to the non-standard terms, we include the leading standard terms. 














































0(s/4M^)°. We assume that the Higgs boson mass is small against the scattering energy, 

Since the standard contributions to the amplitudes W1++++ and _are very large 

(they diverge as cos 6* —> —1), we also include the terms which arise when the sum of the 
leading standard contribution and the subleading (i.e. hrst order non-leading) non-standard 
contribution to these amplitudes is squared, i.e. 2Afstandard-^subieadmg- These terms are 
necessary to describe the effects linear in the for the cross-sections utt and dehned 
in (^). The reason for this is that corresponding leading terms are absent. 

If the couplings are larger, a* > 1, also other sub leading terms might give sizeable 
contributions. Of the possible non-standard subleading terms those which are quadrilinear 
in the couplings will be the largest ones. We include also these terms. They only appear in 
amplitudes with an odd number of longitudinally polarized vector bosons. 

The cross-sections for pp ^ WZ WZ and pp —> Wy —> WZ, calculated with 
the exact (numerically evaluated) expressions for the cross sections for WZ ^ WZ and 
ITy —>■ WZ, respectively, do not deviate by more than 16% from the same cross-sections 
calculated with the high-energy approximation for the WV -^WZ cross-sections presented 
here if the WV invariant mass is in the range 0.8 TeV < a/s < 2 TeV. This is true for 
couplings in the range \ai\ < 0.1. This result was obtained for a pseudorapidity cut of 
7] = 1. For Wy production a similar result can be expected. 

We give expressions for the integrated cross-sections summed over the helicities of the 
outgoing particles. We write the cross-sections as 
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The quantities Ghih 2 h 3 h 4 are the squared helicity amplitudes integrated over the scattering 
angle cos^. 


G 


hl/l2fc3^4 
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j \Mhth2h3hi{cos0p dcosO, 


(43) 


-20 


G is a coupling factor which is different for each process and zq is an integration limit for 
I cos0| determined e.g. by a cut. The indices hih 2 h^hi denote the helicities of the particles 







WV 2 —>■ WV 4 (in this order) and M.hih 2 h 3 hi is the scattering-amplitude. We further defined 
the sums over the transverse helicity amplitudes, 
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The expressions for GT , and GT_ are the same for all processes ITV 2 — WV 4 , = 7 , Z, 
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p and q are the magnitudes of the three-momenta of the vector-bosons in the initial and 
in the final state, respectively, evaluated in the center-of-mass system of the vector-bosons, 
given by 


2 V s 
where i = 2 for p and i = 4 for q. In 
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and tiY = sw/civ- The coupling factors G are given by 
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B.l WZ 


WZ 


For the process WZ ^ WZ there are 25 different helicity amplitudes Mhih 2 h 3 h 4 , which 
cannot be related to each other by discrete symmetries (G, P or T). Of these amplitudes, 
15 have leading terms of the order 0(s/4M^)°, 0{ai) or 0{aiaj). Of these 15 terms, 6 
only appear in the sums and G^_. The remaining 9 integrated squared amplitudes are 
given by 
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where we introduced the variable 
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The subleading terms which are quadrilinear in the couplings are given by, 
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B.2 w~f WZ 

For IF 7 WZ the cross-sections axi and an vanish. There are 27 different amplitndes, 
out of which 14 have terms of the leading order. Of these amplitudes, 8 only appear in the 
sums of the transverse amplitudes and G^_. For the remaining 6 helicity combinations, 
the integrated squared amplitudes are given by the expressions. 
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The non-leading terms for G++++ und G+_have the following expressions. 
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The subleading terms which are quadrilinear in the couplings are given by, 
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4(uvv<i> + c^avv)c^y(uvv<i> iw^B^){zo ^ 

;^3 -I 

+4c^(avi/$ — tw'^B<s>y{zo — 


G^o +00 — ^^wc^ra^/^a^r{9zo — 


B.3 WZ 


W-f 


The terms Ghih 2 h 3 h 4 for ITZ —*> IT 7 can be obtained from the terms Ghih 2 h 3 h 4 for hhh —^ 
ITZ by exchanging the helicity indices according to Also a parity 

transformation ~ might have to be applied. The combinations 

with /i 4 = 0 vanish. 
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B.4 W’’) —^ VF 7 

For IF 7 —IF 7 the cross-sections utl and and the terms Ghih 2 h 3 h 4 with = 0 
vanish. There are 12 different helicity combinations which can not be related to each 
other by discrete symmetries. 8 combinations receive leading contributions. 6 of the latter 
combinations only enter the expressions G++ and G^_. The remaining two leading terms 
are given by 


1 ^ 0 + 0 + — 

Go + o - = + o-b^T{zo + -Zq). 


(56) 


The subleading terms for G 

G-\ —^—I—subleading 


and G_|_are given by 


HwCw 16(avi/<i. aB^)i4:Zofi - lui) -h 16avi/(lni -zq) 


-|-16atyT ®b<i>)(2^o — li^i) — — li^i) T 8(avi/<j) -|- In^ 

G_| - h, subleading = 4c(^/iVl/ + ^5$) 16(2:o/l — lUi) -|- 17Zq + . (57) 

The subleading terms which are quadrilinear in the couplings are given by, 

1 


G++0+ 

0 +++ 

G++0- 

0 +-- 

G+-0+ 

0 ++- 

G+.o- 
0 +- + 


+ o-b^ + o-w) (zo — -^Zq] 

0 




4a^(city<j>-|-ci^<j>)(^0 3^o) 


+4a^(avi/# + o,B^)‘^{zo — -Zq) 


(58) 
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